Resonance states in 12 C and a-particle condensation 
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The states with J n = + , 2 + , and 4 + of 12 C with excitation energies less than about 15 MeV 
are investigated with the alpha condensate wave function with spatial deformation and by using the 
method of ACCC (analytic continuation in the coupling constant) which is necessary for a proper 
treatment of resonance states. The calculated energy and width of the recently observed state 
are found to be well reproduced. The obtained 2^" wave function has a large overlap with a single 
condensate wave function of 3a gas-like structure. The density distribution is shown to be almost 
the same as that of the Oj state that is regarded as a 3a Bose-condensed state, if the energy of 
the 2^" state is scaled down to the same value as the one of the 0^" state. Furthermore, the kinetic 
energy, nuclear interaction energy, and Coulomb interaction energy of the calculated 2^~ state are 
shown to be very similar to those of the 0^~ state. We conclude that the 2j state has a structure 
similar to the 0^" state of Bose-condensate character with a dilute 3a gas-like structure. In addition 
the resonance states, 0^", O4 , 4J, are also discussed. 



I. INTRODUCTION 

In our previous paper we made the strong sug- 
gestion that 12 C, 16 0, and other self-conjugate An nu- 
clei should have a gas-like structure of alpha particles 
in their excited states near the na breakup threshold, 
and that this structure could be regarded as a kind of 
Bose-Einstein condensate, though only a few Bosons are 
concerned in nuclear systems. The j" state of 12 C, called 
the Hoyle state and observed at 0.38 MeV above the 3a 
threshold, and the Ojj" state of 16 0, observed at 0.44 MeV 
below the 4a threshold, were indentified as candidates of 
the three-alpha condensed state and a possible candidate 
of the four-alpha condensed state, respectively by an in- 
vestigation where the following type of wave function was 
adopted: 

$, iQ oc4 [J]exp(-2^|)^(a i )]. (1) 

i=l 

Here ^(a<) oc exp[-(l/26 2 ) £}(r iW - Xi) 2 } and X, = 

"Y^j r j(i) /4, with = 4(i — 1) +j, are the internal wave 
functions and center-of-mass coordinates of the i-th alpha 
cluster, respectively. This wave function has the charac- 
teristic feature that all center-of-mass motions of the con- 
stituent alpha clusters occupy the same spherical 5-wave 
orbit, exp(— 2X 2 /B 2 ). It also distinguishes itself from 
other a-cluster models by the fact that the alpha clus- 
ters do not have any definite geometric configuration, but 
are allowed to move freely with small inter-alpha correla- 
tions or interactions, forming a dilute gas-like structure. 
It should be noted that the above mentioned character- 
istics can be put forward only as far as the variational 
parameter B takes a sufficiently large value so that the 
action of the antisymmetrizer, A, is weak and almost 
negligible. 

It is to be noted that the excited state of a self- 



conjugate An nucleus which has a well-developed na clus- 
ter structure around the na breakup threshold has almost 
zero inter-alpha binding energy and is only held together 
by the Coulomb barrier. Since a particles are Bosons, the 
gas-like state with almost zero inter-alpha binding energy 
is naturally considered to have Bose-condensed structure 
(see also 0). 

For the n—3 case, i.e. 12 C, detailed analysis was made 
by several authors about a quarter century ago, and they 
revealed the fact that several excited states observed in 
the vicinity of the 3a particle threshold have an appar- 
ent 3a cluster structure [3J, |j, |g, |(j . Especially, complete 
3a calculations with fully microscopic models were per- 
formed by two groups, Uegaki et al. £|, and Kamimura 
et al. The calculated binding energy of the 0^" state 
was not much different from the experimental value which 
is 0.38 MeV above the 3a threshold. The 0^" state was 
shown to be dominated by the channel, 1 = 1 = 0. Here 
/ and I indicate the spin of 8 Be and the orbital angular 
momentum of the relative motion between 8 Be and al- 
pha particle, respectively 0. This is consistent with the 
result of the semi-microscopic calculation with OCM (or- 
thogonality condition model) made by Horiuchi, where 
three alpha particles weekly interact in relative S'-waves 
[g. The Oj" state was shown to have a large reduced 
a-width, 9 2 {a) = 0.68 in units of the Wigner limit, 7^ 
= 2>h /2/ia 2 , when the channel radius, a = 6 fm, for 
1 = 1 = 0. Here /1 is the reduced mass between 8 Be 
and the a particle and 6 2 (a) = ^/j^y, where j 2 is the 
reduced a-width (see Eq. (|29[) of section Hi E|) . 

According to the suggestion made in Ref. [| men- 
tioned above, we investigated the 0^ state of 12 C by us- 
ing a deformed alpha condensate wave function which 
only slightly deviates from the spherical one, Eq. 7]. 
We found that each of the 0^~ wave functions, obtained 
by the full 3a microscopic calculations of Refs. Q and 
j5| , has a large squared overlap value of more than 90 % 
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with the single condensate wave function of 3a gas-like 
structure. This result gives decisive theoretical evidence 
that the 0^" state can be understood as a 3a condensed 
state. 

The present paper will mainly focus on the exploration 
of the excitation modes of the finite alpha condensate 
by studying the states of 12 C with quantum numbers 
= + , 2 + , and 4 + with excitation energies less than 
about 15 MeV. The exploration is made by using an a 
condensate wave function with spatial deformation. In 
spite of the fact that the existence of the 2^~ state has 
been suggested at around 3.0 MeV above the 3a thresh- 
old by the former theoretical works Q- 0| ■ it was only 
quite recently that the 2j state was observed at 2.6 ±0.3 
MeV above the 3a threshold together with the alpha de- 
cay width, 1.0 ± 0.3 MeV 0. Morinaga expected the 
existence of a rotational band with the 0^~ state as the 
band head and assumed a 3a linear chain structure as 
the intrinsic structure of the band Q. However, as al- 
ready mentioned above, this assumption is now hardly 
acceptable because the 3a structure of the 0^" state is 
far from that of a linear chain configuration of 3a parti- 
cles. The 2^~ state, in view of the microscopic 3a cluster 
model calculations Q, is dominated by the channel, 1 = 
and I = 2, and has a large reduced a- width amplitude, 
2 (a) = 0.63, with the channel radius, a — 6.0 fm. This 
fact indicates not only that this state has a well devel- 
oped 3a cluster structure like the 0^" state, but also that 
the 2^~ state will be able to be classified in the same fam- 
ily as the 0^~ state in which the single channel, 1 = 1 = 
is dominant. Besides the 2^~ state, the 0jj~ state was ob- 
served at 3.0 MeV above the 3a threshold together with 
a broad width, 3.0±0.7 MeV [lOj. However, there are no 
calculations that reproduce the energy of this state satis- 
factorily. Though the binding energy is higher by about 
3 MeV than experimental observation, the 0^ state ob- 
tained in Ref. Q has large reduced a-width amplitudes 
of both channels, [I <g> I] = [0 <g> 0] and [2 ® 2]. The den- 
sity distribution of the 0g~ state calculated by the AMD 
method (antisymmetrized molecular dynamics) appears 
to be somewhat similar to a linear chain configuration 
of three alpha particles ^lj . Recently a very similar re- 
sult was reported about the 0^" state by FMD (Fermionic 
molecular dynamics) calculation |12| . 

In order to explore the excitation modes of the 3a con- 
densate, we study the 2^~ and 0$ states by a GCM (gen- 
erator coordinate method) calculation in which deformed 
3a condensate wave functions are superposed. This type 
of GCM calculation was successfully used in our previ- 
ous paper 7] for the study of the 0^" state. However the 
GCM calculation was not successful for the 2^~ and 0^" 
states, because this GCM calculation is a bound state 
approximation in spite of the fact that these states are 
resonances with rather large widths of more than 1 MeV 
while the width of the 0^" state is very small of only a 
few eV. Thus, in order to overcome this defect, we apply, 
beyond the bound state approximation, the method of 
ACCC (analytic continuation in the coupling constant) 



to our previous formalism. This method was first pro- 
posed by Kukulin et al. |l3j, and is now recognized as a 
powerful method to calculate the energy and total width 
of a resonance state 0113. 

The main result of the present paper is that the mi- 
croscopic 3a condensate wave function with spatial defor- 
mation is able to reproduce consistently both energy and 
total width of the recently observed 2^ state by the appli- 
cation of the method of ACCC. Furthermore, the 2^ state 
is found to be mostly represented by a single 3a conden- 
sate wave function with slightly deformed shape, where 
the resonance wave function of the 2^~ state is obtained 
by the use of our newly devised technique which allows us 
to extract correct resonance wave functions even within 
the formalism of bound state approximation. The den- 
sity distribution is shown to be almost the same as that of 
the O2 state which is considered as a 3a Bose-condensed 
state, if the binding energy of the 2\ state is adjusted to 
the same value as the one of the 0^" state. The kinetic 
energy, nuclear interaction energy, and Coulomb interac- 
tion energy of the calculated 2 state are further shown 
to be very similar to those of the 0\ state. These results 
indicate that the 2^ state has a structure closely related 
to the 0^~ state where three alpha particles form a gas- 
like structure with dilute density. On the other hand, we 
found that our present functional space of the a conden- 
sate wave function in which strongly prolate condensed 
wave functions are not included is, presumably, not well 
adapted to represent the Ojj" state. This point of view 
may be reasonable if the 0^" state has 3a linear-chain- 
like structure. 

The content of this paper is as follows. In section ITU 
the formulation of the present aproach is explained. A 
brief description of our GCM approach in which the de- 
formed condensate wave functions are superposed after 
angular momentum projection and the complete elimina- 
tion of the center-of-mass spurious components is given. 
The ACCC method is also outlined there. The results of 
the GCM calculations for the 0^~, 2^, 4^, and 0^" states 
are given in subsection 1111 Al The states are treated 
within the bound state approximation. In subsection 
MI Bl the calculations of the energies and widths of the 
2\ and states are presented in which the method of 
ACCC is used. The wave functions of the 2^ and 4^" 
are analyzed in the following subsection IIII CI There we 
introduce a new approximate procedure to estimate the 
wave functions of the resonances. In subsection IIII Dl 
we calculate the energies and widths of the 0^" and O4" 
states by using the method of ACCC. We make the same 
analyses for the wave functions as in subsection IIII Cl In 
subsection IIV Al we check how much = + and 2 + 
components are contained in the deformed intrinsic wave 
function of the alpha condensate. In subsection IIVBI 
we show that if only the orthogonality between the 0^~ 
and 0j states is imposed, the structure of the 0^" state is 
not influenced by the detailed structure of the 0+ state. 
In subsection III CI we demonstrate that the oblate and 
prolate deformed intrinsic wave functions of the alpha 
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condensate give similar wave functions after the angular 
momentum projection. Conclusions are given in section 
El In the APPENDIX, we give the explicit calculus of 
the density distribution in which the center-of-mass spu- 
rious components are completely eliminated. 



II. FORMULATION 

A. Deformed intrinsic not condensed wave function 

The deformed na condensed wave function was first 
proposed in Ref. as a natural modification of the 

spherical shape of wave function, Eq. In this sub- 

section, we explain how to derive the angular momentum 
projected wave function from the deformed na condensed 
wave function. The deformed na condensed state can be 
expressed as given below, taking as a creation opera- 
tor of an a particle with tri-axial variational parameters, 

P x ! Py , Pz ) 



x$ B (# l7 - ••,#„), 



(6) 
(7) 



It is easy to perform the above integration over the vari- 
ables, ili, • • • , Rn, as shown in Ref. ^(| in detail. Con- 
sequently the deformed na-condensed wave function is 
simply given by 



n 

$ na (p) oca [n cx p{/ 

i=l 



2Xf 



( IP 



2X1 



\<I>{oh) 



Here, B 2 



1$ 



(2) 



where 



x(p 0s (n - H)aJ. lTl (ri) • ■■tp 0s (r4 - H)4 4T4 (r 4 ). (3) 

In the above equation, a\ T (r) is a creation operator of 
one nucleon with spin a, isospin r at a spatial point r, 
and (pos (r — R) is a Os harmonic oscillator wave function 
around a center R with the size parameter b as follows, 



(8) 

b 2 + 2p 2 , (k = x,y,z). The other notations 
are the same as used in Eq. Q • A is the antisymmetrizer 
operating on all nucleons in the system. It should be 
noted that Eq. JHJ shows that the center-of-mass mo- 
tions of na particles occupy the same deformed orbit, 
exp{-2X 2 /Bl - 2X 2 /B 2 y - 2X 2 /B 2 Z ). It is also to be 
noted that the normalized na condensed wave function, 
Eq. (JSJ coincides with the shell model Slater determinant 
when (5i —* (i = x, y, z) as a limit. This is an important 
property that is taken over from Brink's wave function. 
On the contrary, the na condensed wave function corre- 
sponds to a free na particle state in which na particles do 
not correlate with each other when /3j — * oo (i = x, y, z). 

The a condensed wave function with good quantum 
number of angular momentum is derived by projecting 
out the angular momentum from the deformed a con- 
densed wave function in the following way, 

$i a (p) = [ dnD(* K (n)R(n)^> na (0) 



<Po. 



(r-R) = (nb 2 r^e, V {- {r 2b f}. (4) = / / • ■ ■ X 



In order to write down the wave function of the de- 
formed na condensed state in coordinate space and 
also to perform numerical calculations of matrix elements 
of physical quantities, it is neces sary to introduce Brink's 
na cluster model wave function [13, <j> B (I?i, • • • , R n ) de- 
fined as 



6X P(-E E ^)* B (^ 1 (^) J Rl,---,i?- 1 (^n) 

i— 1 k—x,y,z ™ 



exp 



$ B (.Ri, •••,#„) 



1 



{ 



det{<A) S (T-i - Ri)Xo 



E E 

2—1 k—x,y,z 



Pi 



•■,21,0, 



VP 



(9) 



(5) 



Here Xar is a spin-isospin wave function of a nucleon. 
From Eqs. 10), ©, the deformed na condensed wave 
function can be written as a superposition of the Brink's 
na cluster model wave functions multiplied by a single 
deformed Gaussian wave packet of R\ , • ■ ■ , R„ with the 
width, /3j, (i = x, y, z), in the following way, 

®na(P) = (riaiTi,---r in a4 n Tin\^na) 

r r " B 2 R 2 R 2 , 



where is the Euler angle, D J MK (Vt) the Wigner D- 
function, R(fl) the rotation operator, and the 3x3 

rotation matrix corresponding to R(fl). Here use is made 
of the relation (po s (R(Q)r — R) = (fos(f — R~ 1 (fl)R). 
We should note that the rotation with respect to the to- 
tal angular momentum is equivalent to the rotation with 
respect to the total orbital angular momentum, since the 
intrinsic spins of the a clusters are saturated. 

Throughout the present paper, all calculations are 
performed with the restricton to axially symmetric de- 
formation in the same way as our previous work, Ref. 
|16| . where the z-axis is taken as the symmetry axis: 



Px = Py 7^ @z- Thus, the formulas used in Eq. © 
are simplified as follows, 



Jdfi 
D J M * K (Cl) 

R(SY) 



Ait J dcost 

Ry(9) 
Ry{6). 



(10) 



In all the following sections, we make use of the notation, 
f3 to express (f3 x = /3 y ,(3 z ). 



B. Elimination of spurious center-of-mass motion 

The total center-of-mass motion can be easily sepa- 
rated out of the na condensed wave function JHJ as fol- 
lows, 



$„ Q (/3)«exp{ - (X Gx + X Gy ) - X G M n M 

X z 

(11) 



2n 



k—x,y,z "* 



(12) 



Clearly Q na ((3) does not contain the total center-of-mass 
coordinate, X G , of the system and it is the internal wave 
function of $ nQ (/3). All the calculations are not made 



with 3>„q.(/3) but with $ na (/3) which is an eigenstate of 
total momentum with eigenvalue zero. Thus, the proper 
a-condensed wave function with a good quantum number 
of angular momentum is not (/3) represented by Eq. 
© but $i a (/3) defined as 



*L(P) = J d cos 9 d J M0 (9)R v (9) $„ 



,09). 



(13) 



In later subsections we use the notation ^J^j' 7 (/3) in order 
to express the normalized state of <I>3 Q (/3) 



09) = 



1> 3 J o(/3) 



<^09)|*LG9)) 



(14) 



Since in Ref. |l6j |. we have given a detailed explanation 
of how to derive the matrix elements of a general transla- 
tionally invariant scalar operator O using the wave func- 
tion with good quantum number of angular momentum, 
$^ Q (/3), we only give a brief outline here. While in Ref. 
jl6j, we examined only the diagonal components of the 
matrix elements in the parameter space (3, we, in this pa- 
per, need to refer to the case of the off-diagonal compo- 
nents. We can write down the matrix elements in which 
the contribution of the spurious center-of-mass motion is 
eliminated, as follows, 



J 



fe(/3)|Q|i>;{ a (/3')) 

\l W a (-9) I Wa 09)) < H a W) I H a m ) 

= J d cos 9 dj (9) ( $ na ((3) | OR y (9) | $ na ((3')) 

y/ J dcos9di (9)($ na (p)\R y (9)\$ na ((3)) ■ J dcos9dJ Q {8)($ na (l3 , )\R y (9)\$ na (f3')) 

= J d cos 9 d£ Q {9){^> na (J3) | ORy (9) | $» a (J3')) [Pp (g) 

sj J dcos9d^(9)(^ na (f3)\R y (9)\^ na {l3))/P 1 {9) ■ J dcos9 d^{9){^ na {l3')\R y {8)\^ na {(3')) / P 2 {8) 

(15) 



Here the functions Pq(9), Pi(9), and P2(9) stem from 
the degrees of freedom of the center-of-mass motion and 
can be easily calculated analytically. We give the explicit 
form of Pq (9) below. The expression of Pi (9) is given by 
simply changing B' x , B', and B' z that appear in the fol- 



lowing expression of Po(9), to B x , B y , and B z . And the 
expression of P2 (9) is similarly given by simply changing 
B x , By, and B z that appear in the following expression 
of P o (0), to B' x , B' and B' z , 



,(9) = (exp(- 
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k—x,y,z 



2n 



Gk 



R y {9) 



exp 
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2n 



B 
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B 2 B' 2 B 2 B' 2 B 2 B' 2 

■^x^x y y z z 



{B 2 y + B' 2 ){{B 2 X + B' 2 )(B 2 + B' 2 ) + {Bl - B 2 )(B' 2 - B' 2 ) sin^ 9} 



(16) 
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C. Approximate equivalence between prolate and 
oblate deformed wave functions after angular 
momentum projection with the exception of 
strongly prolate deformation 

The na condensate wave function with good angu- 
lar momentum $^ Q (/3) has a very characteristic prop- 
erty that <£>^ Q (/3) obtained from a prolate intrinsic state 
can be obtained approximately from an oblate intrinsic 
state and vice versa, except for the case of strongly pro- 
late deformation. In FIG. we show the squared over- 
lap, |($£ J=0 (/3)|$3 a J= °(/3i))| 2 , m the parameter space, 
(3, where /3 X = (5.7 fm, 1.3 fm). FIG. □ shows that 
the + wave functions with prolate and oblate deforma- 
tion are nearly equivalent to each other. On the con- 
trary to FIG. d we show in FIG. [^J the squared overlap 
\(^ a J=0 (j3)\^ a J=0 (/3 L ))\ 2 in the parameter space of (3, 
where /3 L =(0.1 fm, 4.0 fm) which expresses a strongly 
prolate deformation. This figure shows that the angular- 
momentum-projected wave function obtained from a 
strongly prolate deformation with (3 x =(3 v <0.b fm can 
never be similar to any projected wave function from 
oblate deformation. 

We show below that the seemingly strange behavior of 
FIG. n is a generic results in na systems. 



where the Jacobi coordinates, ^ and Xq (i = 1, • • • , n- 
1), satisfy the following relation, 



1 v^ n Y 
n 2~,j = l A 3 

Th (i = !,•••, n- 1). 



Furthermore, in Eq. (|T7Jl . the exponent, (£f x + ) /B 2 - 




£ 2 Z /B 2 , can be rewritten as 



\3B 2 3B?/ Sl \W 2 W 2 



where £j is the polar angle of Thus, the exponential 
term in Eq. I|17|) . i- e - the relative coordinate parts of the 
wave function between n alpha clusters, can be expanded 
as follows: 



Squared overlap (/=0 + ) 



Squared overlap (/=0 4 




4 6 
P x =P y [fm] 




4 6 
P x =P y [fm] 



FIG. 1: Contour map of the squared overlap, 
|{$^ J=0 (/3)|$^ J=0 (/3i)}| 2 , in the two-parameter space, 
@ x {= f3 y ) and p z , where p 1= (p x =/3 y =5.7 fm, /3 Z =1.3 fm). 
The harmonic oscillator size parameter fe=1.35 fm. 



At first, we rewrite Eq. (|12fl in terms of the Jacobi 
coordinates, 

n-i z-2 , <-2 p i 
Siz) 



8 na (/3>^4 exp|-2^/i ! 

i=l 



e 2 +e 



B 2 



bV 



(17) 



FIG. 2: Contour map of the squared overlap, 

I(^3 £ ; J=0 (/ 3 )|S3 Q ' 7=0 (/ 3 l)>| 2 > in the two-parameter space, 
/3 X (= f3 y ) and p z , where p x ={f3 x =f3 y =0.1 fm, /3 Z =4.0 fm). 
The harmonic oscillator size parameter 6=1.35 fm. 
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V 352 
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! = 1 



(18) 



In the above expansion, the parameter, (3B' 2 )^ 1 — 
(SB 2 ) -1 , takes the opposite sign between prolate and 
oblate deformation so that the overlap between the wave 
functions, Eq. i|17|) with prolate and oblate deformation, 
cannot become large. However, once the wave function 
with good angular momentum is projected out of the 
wave function of Eq. (|17|) . the above statement is no 
longer true. From the expansion of Eq. (|18fl . the angular 
momentum projected wave function of the + state takes 
the following form for the leading term, 



A 



exp 



35? 



1 



n—x 

E>#}a< 

»=i 



(19) 



Also for the 2 + state, the leading term is 



A 



exp 



2 1 n ~ 1 

n-l 



X 



£>€?yao(W(a) ■ (2°) 



Eqs. i|19|) an d l|20|l show that the angular momentum- 
projected wave function of the + state as well as the 
2 + state can be parametrized by the sole coefficient 
2(3fl 2 .)-i + (Wl)- 1 . Thus, it follows that the angular 
momentum-projected wave functions from prolate and 
oblate intrinsic states cannnot become so much different 
from one another but rather they are close to each other 
as long as the values of 2(3B 2 .) -1 _|_ (3.BJ) -1 stay similar. 
It is to be noted that this analysis is true for a general 
na system in as much as the leading terms of Eqs. (|19H . 
(|2"U|) are dominant. 



The solution of the above equation leads to the wave 
function as 



(22) 



In this paper, we do not superpose prolate deformed 
wave functions but only §i a ((3) whose /3's satisfy the re- 
lation 0.5 fm <(3 Z < (3 X = (3 y < 20.25 fm. The numbers 
of the adopted wave functions are 61 and 75 for J n = + 
states and the other J w states, respectively. The reason 
of this choice of our functional space is that, as discussed 
in section III CI the prolate and oblate deformed wave 
functions are similar except for the case of strongly pro- 
late deformation. The superposition of the similar wave 
functions sometimes causes numerical inaccuracy. We 
can confirm that the restricted functional space which we 
adopt is sufficient to reproduce the 0+, 2+, and 0^" states. 
For these three states, we solved Eq. I)21f) in our previous 
paper of Ref. , where both the oblate and prolate wave 
functions including strongly prolate ones are superposed. 
We will see that the binding energies and r.m.s. radii 
of these states obtained in the previous paper are com- 
pletely the same as the present results given in TABLE 
III of section IHTAl 

For the 2^ state, we have found that the density dis- 
tribution and the transition density between the 2^~ and 
0^~ states obtained in this paper are almost equivalent to 
those calculated by the wave functions obtained by solv- 
ing microscopic full-three-a problems without any model 
assumptions |5|. We will give the detailed discussion 
about this finding in our next paper. 



E. Analytic continuation in the coupling constant 



D. Hill- Wheeler equation 

In this subsection, we give a brief outline of the GCM 
approach. In order to discuss the quasi bound states as 
well as bound states, we solved the following Hill- Wheeler 
equation by superposing the wave functions of Eq. I|13|) 
with good quantum numbers of angular momentum in 
which the center-of-mass spurious components are com- 
pletely eliminated, 

]T ($L(P)\(h - m\*LW)tiW = o- (21) 



Besides the 0^" state which has a narrow width, other 
excited states located above the 3a threshold have non- 
negligible widths. In order to discuss resonance states 
that have broad widths, we go beyond the framework of 
bound state approximation by using the ACCC method. 

In the ACCC method, we introduce an attractive 
pscudo potential, V, which is added to the microscopic 
Hamiltonian H in the following way: 

H'(6)=H + SxV, (23) 

where S is a newly introduced coupling constant. The ba- 
sic idea of the ACCC method is to distinguish resonance 
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states from continuum states by increasing the coupling 
constant 5 from the physical value, 5 = 0, and making 
the interaction unphysically attractive. As a result of this 
operation, the energy of a resonance state goes below the 
3a threshold, while a continuum state never goes under 
the 3a threshold. This operation allows us to handle the 
resonance states within the framework of the bound state 
approximation. 

In this work, we adopted the following form for the 
pseudo potential, 



V = 



(24) 



The reason why we assign 2.5 fm for the range parameter 
of the pseudo potential is that in order to pull down the 
resonance state efficiently into the bound state region, 
it is useful to make the range parameter of the pseudo 
potential larger than the range of the attractive part of 
the nuclear force. 

In order to treat the resonance states which are con- 
verted to bound states, we solve the following Hill- 
Wheeler equation for the ACCC Hamiltonian H'(8), 

J2 ($L(P)\(H'(5) - E{(6)) |$ 3 J Q (/3')) //(£', 8) = 0. 

(25) 

The solution of the above equation is given as a function 
of the coupling constant 5 as follows, 



*a / W = E/a(/3,'5)$3 J q (/3)- 



(26) 



Once the resonance state is obtained in the bound state 
region, the energy as a function of 5 of the state in ques- 
tion can be approximated by a Pade approximant that 
has the following form, 



k 



[N,M] 



. C + C\X + C2X + . . . + CnX 



N 



E(S) - El h a (5) = k? 



1 + d\x + d,2X 2 
2 

[N,M] ' 



d M x 



M ' 



(27) 
(28) 



Here x = ^5 — So, and the constant So satisfies the rela- 
tion, E(S ) = E^(S ) = 3E a (S ). E a {S) is the binding 
energy of an alpha particle as a function of 5. We can 
determine the coefficients of the Pade approximant of the 
right hand side in Eq. l(27|) using the values of the bind- 
ing energy derived from the Hill- Wheeler equation, Eq. 
(|25|l . While the energy function, k 2 [N, M], determined 
with the help of the Pade approximation only has a real 
part in the bound state region below the 3a threshold, 
it takes complex numbers above the 3a threshold, where 
the binding energy and the total width will be given as 
Re{E(S) - E&(6)) and -2 • ha{E(S) - E^(S)), respec- 
tively. We can extrapolate the energy function in <5 back 
to the physical point, S = 0, to reach the desired binding 
energy and width. 



The numerical accuracy of the binding energy which 
we obtain by the GCM calculation is approximately seven 
significant digits, namely the precision of several tens of 
eV since the binding energy is around 80 MeV. The deter- 
mination of the coefficients Cj and di of the Pade approxi- 
mant is made by using E(S) — E^ which is of the order of 
10 _1 MeV ~ 1 MeV. Therefore, due to this subtraction 
of E(S) — E^, the numerical accuracy of the coefficients 
Ci and di is lower than that of the GCM binding energy 
and it atteins approximately four significant digits. If the 
width of a resonance is much smaller than the energy of 
the resonance by more than four orders of magnitude as 
this is the case for the state whose energy and width 
are about 0.4 MeV and about 9xl0~ 6 MeV, respectively, 
we cannot calculate the width of the resonance with the 
ACCC method because of lack of numerical accuracy. 

However, we can devise an approximate method to cal- 
culate such small widths within the framework of the 
ACCC method by using a region of negative values of 
S where we can safely calculate the resonance width 
T(S). With the condition that the width be calculated 
with good numerical accuracy, we choose a negative S, 
S=5 Z <0, whose magnitude |<5 Z | is as small as possible. 
According to the i?-matrix theory 0, the decay width 
r(<5) is expressed as 

T(5) = 2P l (e(S))- 1 2 (8), 

Pl{e) = F?( P ) P +Gi( P y p = a \f¥ e > (29) 

where P; is the penetrability of the 8 Be(0 + )+a decay 
channel with partial wave /1 being the reduced mass 
of 8 Be and a, and a the channel radius. Here e(S) is the 
decay energy, e(5)=E(S)— E^=° (S) , where E^°(S) is the 
binding energy of 8 Be(0 + ), and Fi and Gi are regular 
and irregular Coulomb functions, respectively. From this 
formula we obtain 



T(S = 0) 



T(S Z 



(30) 



under the approximation 7 2 (<5 = 0)~7 2 (<5 2 ). Here, e is 
the observed decay energy into the 8 Be(0 + )+a decay 
channel. 



F. Microscopic Hamiltonian 



The microscopic Hamiltonian is expressed as 



H = T -T G + V N + Vc. 



(31) 



Here T denotes the kinetic energy, Tg the total center- 
of-mass kinetic energy, and Vc the Coulomb interaction 
between protons. These quantities are given by 



T 
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(32) 



Vn denotes the effective two-body nuciear interaction. 
In this paper, the Volkov No. 1 force with the value of 
Majorana parameter M — 0.575 used by Uegaki et 
al. and the Volkov No. 2 force with M = 0.59 [3 used 
by Kamimura et al. are adopted, which are expressed in 
the following formula: 



1 12 

1 m 



M) - MP a P T } l0 J2 vn exp 




(33) 



III. RESULTS 



A. GCM study of the ground band states and the 
O2 state 

In our previous paper of Rcf. jj], we concluded that 
the 0^" state of 12 C has Bose-condensate character where 
the three alpha particles occupy the same center of mass 
5-state. Here we give some more detailed analyses to 
reinforce this conclusion. 

In Ref. [7J, we analysed the microscopic 3a calcula- 
tions by Uegaki et al. and Kamimura et al. which were 
performed about a quarter century ago. We then con- 
cluded that the 0^" wave functions obtained by them are 
almost the same as the 3a Bose condensed wave func- 
tions. We here briefly recapitulate the essence of our 
argument in Ref. 0. 

In Ref. 0, Uegaki et al. adopted the Volkov No. 1 
force with a Majorana parameter, M = 0.575, as effective 
nuclear force, while in Ref. p| , Kamimura and co-workers 
adopted Volkov No. 2 force with the value, M = 0.59. 
Uegaki. et al. and Kamimura and co-workers followed 
different approaches, but in essence there was no differ- 
ence between them since both approaches fully took into 
account the degrees of freedom of relative motions be- 
tween the constituent alpha particles as well as the Pauli 
principle. In both approaches, the alpha particle with 
(0s) 4 of harmonic oscillator type was used with a fixed 
oscillator size parameter value, b = 1.41 fm for Volkov 
No. 1 and b = 1.35 fm for Volkov No. 2. For example, 
Kamimura and co-workers adopted the resonating group 
method (RGM) where the following equation was solved, 



(cb 3 (a)\(H-E)\A{ X (s,t)cj ) 3 (a)})=0. 



(34) 



Here s and t are the Jacobi coordinates of the center- 
of-mass motion of three alpha clusters. We can clearly 
see that our model wave function, Eq. (|12fl . mathemati- 
cally corresponds to a specific form of the above complete 
three-body wave function A{x(s, t)<fi 3 (a)}. 
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FIG. 3: Energy spectra of low lying positive parity states of 
12 C. "Full 3a" for the Volkov force, No. 1 and No. 2, are 
taken from Q and 0, respectively. 

(a) : the 2% state in experiment comes from Ref. 0- The 
other observed states are taken from Ref. [Tcj . 

(b) : GCM + ACCC 



Adopting the same nuclear forces as Uegaki et al. and 
Kamimura et al., we solved the Hill- Wheeler equation of 
Eq. ij?l"1l for the ground state band with 0^, 27^, 4+, and 
the O2 state. The eigen energies and r.m.s radii of our 
GCM calculations are compared with those by Uegaki 
et al. and Kamimura et al., respectively, in TABLE [I] 
The comparison of the eigen energies is also shown in 
FIG. Throughout the 0^,2^,47;, and 0^ states, as 
well as our binding energies are almost the same as those 
obtained by Uegaki et at, and also by Kamimura and 
co-workers. Considering the mini-max theorem of the 
variational problem and the fact that the complete 3a 
model space contains our model space, we can say that 
our wave functions of those states are almost identical 
with the ones given by Uegaki et al. and by Kamimura 
and co-workers for each case of the adopted nuclear force. 

The calculated large root mean square radius of the 0^ 
state is quite remarkable. It is 4.44 fm for Volkov No. 1 
and 3.83 fm for Volkov No. 2. The root mean square radii 
of the three states of the ground band are very similar 
and show that these states have a compact structure with 
normal density and do not have a pronounced three-alpha 
cluster structure. 

The most important result in Ref. [7J is the fact that 
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TABLE I: Comparison of the GCM calculation solving the Hill- Wheeler equation, Eq. 12 H . with the full 3a calculation 0,0 
with respect to the binding energy, E, and r.m.s radius, -Rr.m.s- The comparison is made for the 0^", 2f , 4^, and 0^ states 
for two cases of the effective two-nucleon force. The r.m.s radii of Ref. Q which was cited in our previous paper Q were the 
charge radii containing the proton size -J (r 2 ) p — 0.813 fm. Here we cite the values without proton size effect, which can be 
directly compared with our r.m.s radii. Note that for the 4+ state, the spurious components of continuum states are subtracted 
by the method discussed in section fill CI 

Volkov No.l M = 0.575, Ef a = -81.01 MeV Volkov No.2 M = 0.59, Ef a = -82.04 MeV 



GCM calculation full 3a calculation [4] GCM calculation full 3a calculation [5] 



r 


E (MeV) 


-Rr.m.s (fm) 


E (MeV) 


-Rr.m.s (fm) 


E (MeV) 


-Rr.m.s (fm) 


E (MeV) 


Rr.m.s (fl 


of 


-87.81 


2.40 


-87.9 


2.40 


-89.52 


2.40 


-89.4 


2.40 


21 


-85.34 


2.38 


-85.7 


2.36 


-86.71 


2.38 


-86.7 


2.38 


4f 


-79.87 


2.30 


-80.4 


2.29 


-80.30 


2.31 


-80.0 


2.31 


0+ 


-79.97 


4.44 


-79.3 


3.40 


-81.79 


3.83 


-81.7 


3.47 



the wave function calculated in our GCM approach, 
therefore the 0^" wave function given by Uegaki et al. or 
by Kamimura et al., is almost completely equivalent to 
the much simpler condensate wave function. In getting 
this result, the following steps were taken. At first we 
calculate the energy surface after the angular momen- 
tum projection to J 71 " = + and obtained the minimum 
energy point at /3 = (1-5 fm, 1.5 fm) where the minimum 
energy is calculated as —87.68 MeV for the case of Volkov 
No. 2 force. The minimum energy state is considered to 
correspond to the ground state. We then calculated the 
projection operator, 

Pf=°(/3) = 1 - |^/ =0 (/3))(^ J=0 (/3)|, 05) 

at /3 = (3q. Here ^^ J (/3) is the normalized state 
of $3 Q (/3) as defined in Eq. 1(13)1. The operator 
Pf =0 (j3o) creates the functional space orthogonal to 
the approximate ground state $3=°(/3 ). Next we 
calculated the energy surface within this functional 
space spanned by Pf_ =0 (/3 ) < E ) ;^ (/3). The minimum 
energy point of this energy surface was obtained at 
j3 = j3i = (5.7 fm, 1.3 fm) in the case of Volkov 
No. 2 force. Finally we calculated the squared over- 
lap |(^|CiPi =0 (/3 )S^=°(/3i))| 2 withCi standing for 
the normalization constant of P/ =0 (/3 )<i>^ (/3i) and 
we found that this quantity is as large as 0.97. The 
fact that the 0^" state is represented by a simple state 
Pj 7=n (/3 )(f>^= (/3i) naturally leads to the correctness of 
our picture that the 0^" state is considered to be a Bose 
condensed state composed of 3a particles. 

On the other hand, we here analyse the wave function 
of the O2 state, by simply calculating the squared 

overlap, |($ 3 ^ [ J_0 (/3)| V I' A '=2)I 2 ' m tnc parameter space, (3. 
We show the contour map in FIG. 0] These quantities 
indicate how much the Of wave function can be repre- 
sented by a single condensed wave function with good 
angular momentum + . The maximum value amounts 
to 0.82 at (3 = (32 = (6.5 fm, 1.5 fm) for the case of 
Volkov No. 2 force. The same calculation gives 0.88 at 



/3 = I3 3 = (7.5 fm, 1.5 fm) for Volkov No. 1 force. This 
should be compared with the fact that the maximum 
squared overlap between the Of wave function obtained 
by Uegaki et al. and the Brink's + wave function 
with a single three-alpha configuration was at most 0.5. 
This result gives further support for the interpretation 
that the 0^ state has a gas-like structure where three- 
alpha particles are condensed into an indentical iS-wave 
orbit. 
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FIG. 4: Contour map of the squared overlap, 

K$3a J=0 (/3)l*A=a)| 2 j in the two-parameter space, X (= (3 y ) 
and f3 z . Numbers attached to the contour lines are the 
squared overlap values. The adopted effective nucleon force 
is Volkov No. 2. 

We should keep in mind that the wave function 
$ 3 7 t ^°(/3) of Eq. I|13|) necessarily contains compact com- 
ponents like the shell model state even if the size param- 
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eter (3 is large. Therefore the overlap ($^ J ~° (/3) | ^a=°) 
even for large f3 = f3 2 or (3 = /3 3 suffers from the ef- 
fect of the compact components contained in ^i^ 7-0 (/3) 
with (3 — (3-2 or (3%. In order to calculate the overlap 
of w hh the gas-like state which does not contain 

the compact components, we need to subtract the com- 
pact components from $3^' /_0 (/3) . The wave function 
Pf =0 (f3o)^>^°(f3) is just such a kind of wave function in 

which compact components represented by ®ia°{flo) are 
subtracted. As we mentioned above, the squared overlap 
of ^x=2 w ith P/ =0 (/3o) < 1 ) 3q °(/3) can become as large as 
0.97. 

We now discuss the ground rotational band states with 
J n = 0+, 2+, and 4+. As we mentioned already, the 
binding energies of the ground rotational band states by 
our GCM are in good agreement with those by Uegaki 
et al. and Kamimura et al., and therefore our GCM 
wave functions of these states are almost identical to 
their wave functions. It is not unnatural that our GCM 
wave functions could well represent these states, if we 
consider the fact that the normalized wave function of 
Eq. (|12fl goes to a shell model Slater determinant when 
Px = (3 y — (3 Z — * 0. However these states have a quite 
different structure from that of the 0^ state and cannot 
be considered to form a dilute 3a cluster structure. 

In TABLE. [HJ we give the kinetic energy, nucleon- 
nucleon interaction energy, and Coulomb interaction en- 
ergy for each of the Of , 2+, 4^, and 0^ states. We should 
note that the internal contributions of the alpha particles 
are subtracted. Kamimura et al. calculated the kinetic 
energy for each of the 0+, 2^, and 0^ states in Ref. y|, 
and our result well agrees with theirs. We see that for 
the 0^" state, the kinetic and nucleon-nucleon interaction 
energies are much smaller, compared with the case of the 
ground rotational band composed of the 0^, 2+, and 4^ 
states. The Coulomb interaction energy of the 0^ state is 
smaller than those of the ground rotational band states. 
These results are well understood by the fact that the 
0^ state has a strongly developed 3a cluster structure 
with low density where the inter-alpha binding is small. 
Simultaneously, it is indicated that the ground rotational 
band states have a compact structure with normal den- 
sity, as is clearly shown by correspondingly small values 
of r.m.s radii of these states. 

In FIG. |SJ we show the density distribution which is 
the expectation values of density operater defined as, 



p(a) 



a). 



(36) 



A brief explanation of how to derive the expectation 
values of the above density operator will be given in 
the APPENDIX. Obviously the wave function of the 0^ 
state leads to a low density structure with a long tail, 
while the states of the ground rotational band have com- 
pact structures with all almost the same density distri- 
bution. This result provides a strong support to the idea 
that the 0^ state can be described as a dilute density 




FIG. 5: Density distributions of the 0^, 2f, 4f , and 0^" states 
which are the expectation values of the density operator de- 
fined in Eq. (13611 . The adopted effective force is Volkov No. 
2. 



state forming a gas-like structure composed of three- 
alpha particles, as characterized by the wave function, 
Eq. (|12fl . On the contrary, the ground rotational band 
states with J 71 " = 0^, 2^, and 4+ have a quantitatively 
different structure from the one of the 0^" state. Hence 
even though these states could be well represented by 
our GCM wave functions, they could never be described 
in terms of a Bose-condensation composed of three-alpha 
particles. 

The ground band spectrum obtained by us as well as 
those by Uegaki et aL0 and by Kamimura et aZ.JjJ are 
quite compressed in comparison with experiment. This is 
a well-known problem of microscopic cluster model calcu- 
lations. Takigawa and Arima [2(| showed that the effect 
of the spin-orbit force taken into account by broken sym- 
metry states with [4431], [4422], and so on introduced 
into the model space in addition to the cluster model 
states with maximum spatial symmetry [444] can largely 
remedy this shortcoming of the cluster model. The pair- 
ing correlations may also have to be considered. 



B. ACCC study of 2+, 4+ and 

The 2^ state is located at 2.6 ± 0.3 MeV above the 3a 
threshold with the a-decay width 1.0 ± 0.3 MeV. There- 
fore the bound state approximation used for the Oj state, 
which is only 0.38 MeV above the 3a threshold, is no 
more useful for t he 2\ state. This situation is shown in 
FIGs. |6(a)]|6(b) 

we give the contour map of the en- 



In FIG. 6(a) 



ergy surface for the 2 + states in the parameter space, (3. 
We obtained the minimum energy point at (3 = (3^ = 
(1.5 fm, 0.35 fm). The optimum wave function is written 
{(3q), which corresponds to the 2f state, since 



as $3 7 " 2 
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TABLE II: Expectation values of kinetic energy, nucleon-nucleon interaction energy, and Coulomb interaction energy of the 
relative motions between a particles, denoted as T — Tq, Vn, and Vc, respectively. Internal contributions of the a particle are 
subtracted. Calculation is made for the Of, 2f, if, and Of states for two cases of the effective two-nucleon force. Note that 
for the 4+ state, the spurious components of continuum states are subtracted by the method discussed in section llll CI Units 
are in MeV. 
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the binding energy, —84.65 MeV is similar to the correct 
binding energy of the 2f state, —86.70 MeV, as indicated 
in TABLE [I] In order to look for the 2^" state, we define 
a projection operator, 



P[= 2 (f3») = l-\<t> 



u "W?))(^(fl?)l 



3 a 



(37) 



Using this, we calculated the 2 + energy surface of the or- 
thogonalized states, Pf= 2 ((3l?)&£= 2 (f3), to the min imum 
energy state. The contour map is shown in FIG. |6(b) 
There appears no minimum point in this map. This is 
because the 2 J state lies close to the top of the Coulomb 
and centrifugal barriers. This situation of the 2^~ state is 
in marked contrast to that of the Oj" state for which, as 
shown in Ref. Q , we could see a clear minimum point in 
the contour map of the energy surface for the + states 
calculated in this way. This is because the 0^ state lies 
well below the top of the Coulomb barrier. Therefore, 
the 2j" state has to be correctly treated as a resonance 
state. 

In FIGs 7(a)|7(cj we plot the energy eigenvalues of the 
ACCC Hamiltonian H' of Eq. J23|) obtained by solving 
the Hill- Wheeler equation, Eq. fgft, for J n = 0+,2+, 
and 4 + . The solutions of the above equation of motion 
depend on the coupling constant, S, and the energy eigen- 
values are plotted as a function of d. It is seen that res- 
onance solutions dive under the line of the 3a threshold, 
E = 0, as 5 increases. On the other hand, the con- 
tinuum solutions, which are insensitive to variations of 
S and look just like horizontal straight lines on the fig- 
ures, stay above the line of the 3a threshold. This is 
how resonance states can be separated from continuum 
states. It is to be noted that around S — 0, an undesir- 
able mixing of resonance and continuum states occurs, 
since the Hill- Wheeler equation, Eq. (|25|l . has not been 
solved with proper boundary conditions. Therefore the 
correct resonance states cannot be clearly identified. As 
characteristic features of these figures, the trajectory of 
the O2 state as well as the 4+ state can be clearly traced, 
even when the coupling constant, 5, decreases to zero, 
while the trajectories of the other resonance states, for 
instance, that of the 0;j~ state get obscure when 5 becomes 
less than about 0.01. As a result, the 0;j~ state cannot be 



distinguished from a few neighboring states at 6 = 0. It 
follows that it is a very good approximation to solve the 
2 state within the framework of the bound state approx- 
imation as practiced in the previous subsection, while for 



the 2i 



o 3 + , 



and 4 J states with non negligible widths, it 



is necessary to go beyond the framework of the bound 
state approximation. In order to solve those resonance 
states, the method of ACCC is applied here. 



In Figs 8(a) |8(b) we show the converged complex en- 



ergy functions for the 2^ and 4f states. The energy func- 
tions of resonance states have the form of Eqs. I|27l) and 
(|28|l . As already mentioned in section III El once the co- 
efficients of the Pade rational function are determined by 
the use of the calculated energy eigenvalues underneath 
the 3a threshold, the position and total decay width can 
be derived by extrapolating the complex energy functions 
of the coupling constant, 8, to 5 = 0. For the 2jj , and 
42 states, we made use of energy region down to about 
3 MeV starting from several tens of keV below the 3a 
threshold in order to determine the coefficients of the 
Pade approximant. 

In FIGs. |8(a)[ |8(b)| the calculated stable trajectories 
for a few kinds of Pade approximants labeled as [N, M] 
(N = M) defined in Eq. H27J are drawn by superpos- 

7(c)| respectively. We see that 



ing them on FIGs. 7(b) 



the real parts of all calculated trajectories re trace, even 
above the threshold, the trajectories in FIGs. 7(b) 7(c) 
It seems that we need to adopt the Pade approximant of 
N = M > 6 for these states in order to correctly retrace 
the real parts of the complex energy functions, especially 
above the threshold, the energy eigenvalues obtained by 
the bound state approximation. On the other hand, we 
applied the Pade approximations of N = M > 9 to the 
2^~ state, and JV" = M > 8 to the 4 J state, but we did 
not adopt the results of these applications because they 
appear to be ill-behaved or divergent. It should be noted 
that the two kinds of analytically continued Pade approx- 
imants ofdifferent N = M values for the each state on 
FIGs. |8(a)| |8(b)| almost coincide so that we can hardly 
distinguish them from one another, 
state, FIG. _ 
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FIG. 6: (a): Contour map of the energy surface of the 2 + state 
in the two-parameter space, /3 X {= /3 y ) and f3 z . (b): Contour 
map of the energy surface corresponding to the 2 + state or- 
thog onalized to the state at the minimum energy point in Fig. 
|6(a)| In both figures, the adopted effective nucleon force is 
Voikov No. 2, and numbers attached to the contour lines are 
binding energy values given in units of MeV. 
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FIG. 7: The energy eigenvalues of the Hill- Wheeler equation, 
Eq. 125II . plotted as a function of the coupling constant, S. 
(a): For = 0+. (b): For J 7 * = 2+. (c): For J 71 = 4+. The 
adopted effective nucleon force is Voikov No. 2 for all cases. 
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MeV, when the Volkov No. 2 force is adopted. The cal- 
culated binding energy, 1.55 MeV, is quite a bit smaller 
than the observed value, 2.6 ±0.3 MeV. According to the 
lower binding energy which makes the decay penetrability 
smaller, the total decay width is also much smaller, i.e. 
0.09 MeV, than the observed value, 1.0±0.3 MeV. When 
we adopt the Volkov No. 1 force, the energy position is 
calculated as, Re(E(S = 0) - E^(S = 0)) = 2.1 MeV. 
According to the higher energy position, the calculated 
total decay width is larger than the one of the Volkov No. 
2 case and is, r = -2 • hn(E(5 = 0) - E&(6 = 0)) = 0.64 
MeV, which quite well reproduces the experimental data. 
The sensitive dependence of the penetrability on the en- 
ergy position can be seen by adjusting the value of 5 
and by making the whole system a little less bound, so 
that the real part of energy function takes a reasonable 
value. When S = —0.0035, the real part of energy func- 
tion, Re(£ [8 , 8] (<S = -0.0035) - Ef a {5 = -0.0035)) = 2.0 
MeV, and Y = 0.61 MeV. At this value of 6, the thresh- 
old energy of 8 Be + a turns to be 0.172 keV measured 
from the three alpha threshold energy, though at 5 = 
the 8 Be+a threshold energy is —0.244 MeV. When we 
adopt the Volkov No. 1 force, the threshold energy of 
8 Be + a at 6 = is 0.191 MeV measured from the three 
alpha breakup threshold energy. Consequently, we can 
say that for both effective nucleon forces, the resonance 
energy and width of the 2^~ state are mutually consis- 
tent and in good agreement with the recently observed 
experimental data. 

There is another way to adjust the present unrealistic 
penetrability in the case of the Volkov No. 2 force to a 
more reasonable value. It is using the i?-matrix theory 
in which the decay width is expressed as T — 2P\ ■ j 2 , as 
given in Eq. i|29|) . From this relation, we can estimate the 
total decay width corresponding to the observed energy 
position of the l\ state, 2.5±0.3 MeV above the 8 Be + a 
breakup threshold energy which is 0.092 MeV above the 
three alpha breakup threshold. The estimation is made 
through the following relation: 



r' = 



Pi«5 = 0)) 



T(6 = 0). 



(38) 



The calculated total dacay width, V is 0.24 MeV when 
the channel radius a = 6 fm, which is still smaller than 
the observed value. It is to be noted that the above 
expression of penetrability is an approximate one and 
only applicable to the case of a narrow resonance width. 
Since the observed width of the 2^ state which is 1.0±0.3 
MeV, is no longer narrow, the use of this expression of 
penetrability may lead to an unsatisfactory result. 

In the present work, we do not argue about the 2^ and 



2\ states which are shown in Fig. 7(b) A recent exper- 
iment assigns a 2 + state at the excitation energy, ~ 14 
MeV with the width, ~ 1 MeV in the preliminary anal- 
ysis [2lJ. Theoretical predictions that a 2 + state exsists 
in this energy region with the width, ~ 1 MeV are also 
reported |22 
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FIG. 8: (a): Analytically continued complex energy fun ctions 
for the 2 J state are drawn by superposing them to FIG. |7(b)1 
Note that 
tories can 
figure, (b): 
for the 4^ 
Note that 

tories can be hardly distinguished from one another in this 
figure. 
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Analytically continued complex energy fun ctions 
state are drawn by superposing them to FIG.|7(c)| 



k 



[6,6] 



and kyj y\ 



are converged, so that their trajec- 



In section III El we devised an approximate method 
to calculate the width which is much smaller than 
the energy measured from the threshold within the 
framework of ACCC method. We calculated the 
width of the 0^ state by using Eq. iJSDJl, T(8 = 
0)«(P; = o(e)/P; = o(e((5 z )))r((5 z ), where the observed decay 
energy e=0.288 MeV. When <5 z =-0.0035, the calculated 
width r(c5 = 0)«7.0 x 10~ 6 MeV, where the decay energy 
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and width are e(<5 z )=0.643 MeV and r(5 z )=0.012 MeV, 
respectively. When 6 Z =— 0.005, the calculated width 
T(S = 0)«9.9 x 10~ 6 MeV, where the decay energy and 
width are e(<5 2 )=0.790 MeV and r(<5 z )=0.064 MeV, re- 
spectively. The adopted channel radius is a=6 frn in this 
calculation. The obtained values are consistent with the 
observed width 8.7 ± 2.7 eV pjjj. 

In FIG. |8(b)] the energy function of the 4^~ state using 
the Volkov No. 2 force is drawn. The obtained binding 
energy mesured from the threshold and the total decay 
width are 4.5 MeV and 4.6 MeV, respectively. When we 
adopt the Volkov No. 1 force, we obtain 5.4 MeV and 
4.0 MeV for the energy and the width, respectively. In 
the case of the Volkov No. 1 force, Uegaki et al. pre- 
dicted in the Ref. Q the exsistence of the 4^ state that 
is, however, experimentaly unknown. In spite of the fact 
that we obtained a unique and stable trajectory for the 
4^" state for both of the adopted effective nuclear forces, 
we, however, could not make a decisive argument about 
whether the obtained 4^ state possesses a correct phys- 
ical meaning. The reason for this will be clarified in the 
following subsection. 



C. Analysis of wave functions of the and 4^"- 
states 



of the level-level repulsion between the resonance-state 
and continuum-state trajectories is no more existent, we 
may consider that the mixing of the continuum state into 
the resonance state is very small. 
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FIG. 9: Density distributions of the wave functions, 
*(r| (5=0.00542), #^(5=0.00265), *^(<5=0.00165), and 
^x=e(S=0), together with that of the state given in the 
section, 1111 Al The density operator is defined in Eq. 1361 . 
Volkov No. 2 force is adopted. 



We have seen that our wave function of the 2^ state 
reasonably reproduces the observed binding energy and 
alpha decay width, with the help of the ACCC method. 
In this subsection, we analyse the wave function of the 2 J 
state. Since the method of ACCC, unfortunately, does 
not provide the wave function of the resonance state, we 
utilize the GCM wave functions for the ACCC Hamil- 
tonian H'(8) in order to estimate the correct l\ wave 
function. 

FIG. 7(b)| shows that the binding energy of the 2 J state 
coincides with the threshold energy when 5 = 0.00626. 
Let V{(6) be a solution, Eq. (gBJ, of the Hill-Wheeler 
equation, Eq. (|2*5)l . For 6 < 0.0062, as long as we can 
safely trace the clearly visible trajectory of the resonance 



state corresponding to the 2J" state shown in FIG. 7(b) 



we will be able to regard the wave function \I/^ =2 (oj on 
the trajectory of the 2% state as consisting of an almost 
pure resonance component for the ACCC Hamiltonian 
H'(S). Below we give arguments for this conjecture: as 
already mentioned, the horizontal straight lines in Fig. 
7(b) express the trajectories of the continuum states. We 



can check this fact by calculating the density distribu- 
tion p(a) of any state on the continuum-state trajecto- 
ries, which shows that p(a) < 0.05 for a < 10 fm. The 
coupling of the resonance state with the continuum state 
is the strongest at the point where the resonance-state 
trajectory crosses the continuum-state trajectory. The 
coupling of the resonance state with the continuum state 
becomes weaker as we go away from the crossing point 
along the resonance-state trajectory. When we arrive at a 
point on the resonance-state trajectory where the effect 



In FIG. OH we show the density distributions p(a) 
of four states on the 2^~ resonance trajectory, ifa = 



*a 7 =I (^=0.00542), ipjx = *a=3 ((5=0.00265), ip m = 
(£=0 .00165), and ?/>iv = ^\ZK 6 = °)- As is seen in 
FIG. |7(b)j the first two states, tpi and V'n, are out of in- 
fluence of the crossing with continuum trajectories. How- 
ever, the last state, ^>rv> seems to contain non-negligible 
amount of the continuum state component, because this 
state is located near the crossing point between the 2^" 
trajectory and the fourth lowest continuum trajectory 
with large mutual repulsion. In addition, it may be said 
that the state, ipiu, has little mixing with the second and 
third continuum states, since the state is closely located 
between the second and third crossing points. For the 
sake of comparison, we also show in FIG. [21 the density 
distribution of the 0^" state which we already showed in 
FIG. El We note that the density distributions of the first 
two states, tpi and "0n, are confined inside a finite spatial 
region. Actually the integration of p(a) 



I(R) 



p(a)da 



(39) 



which should be unity for R = oo, I(R = oo) = 1, be- 
comes almost unity already for finite R except for the last 
two states, ipm and tpiy. Such finite values of R are 8 fm 
for the 0^" state and and 11 fm for tpu. In the cases 
of tjjni and ipiy, J(i?=13 fm)=0.96 and 0.90, respectively, 
and p(a) is not damped away even beyond a = 15 fm. 

The reason for the choice of S = 0.00542 of the first 
state -01 = ^a=2(^ = 0.00542) is to make its excitation 
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energy measured from the 3a threshold, E(6) — E^(6), 
be equal to that of the 0^" state which is 0.26 MeV for 
our present choice of Volkov No. 2 force. We see that 
the density distribution of this state, ipi, is almost the 
same as that of the Oj" state. This result implies that if 
we pull down the excitation energy of the 2^ state to the 
same value as that of the state, the density distribu- 
tion of the 2^ state becomes almost identical to that of 



the 0J state. In FIG. 10(a) we show the squared overlap 



between tp\ and the single 2 + condensate wave function 
(/3M)| 2 



*£ J=2 (/3) 



- 3a 



($3 a J - 2 (/3M)| 2 . We see that the max- 



imum value of the squared overlap is more than 90 %. 
Here we should recall that in the case of the 0^" state the 
maximum value of such squared overlap with single + 
condensate wave function S^" 7-0 (/3) amounts to 82 % 
for the Volkov No. 2 force and to 88 % for the Volkov 
No. 1 force. We further calculated the squared overlap of 
-01 with -Pj 7=2 (/3o') ( I ) 3 Q =2 which is the projection operator 
defined in Eq. I|37fl . We found the maximum value to be 
more than 95 %. In the case of the 0^ state, a similar cal- 
culation of the squared overlap with Pj 7=0 (/3 )<i>;{^ (/3 1 ) 
gives a value as large as 97 %, where /3 1 =(5.7 fm, 1.5 fm), 
as mentioned in section ITlI Al We thus can say that the 
wave function ipi has a very similar structure as the 0^" 
state of Bose-condensate character. 

We see in FIG. that the density distribution p(a) 
of the second state, ifju=^ J x =l{5=Q.002Qb), is slightly 
pushed outwards compared with that of ip%, which we 
consider to be natural because the excitation energy of 
tpu measuredd from the 3a threshold, E(5) — E^(S), is 
1.01 MeV and hence is higher than that of ipj which is 
0.26 MeV. We also calculated the squared overlap be- 
tween tpu and the single 2 + condensate wave function 

®3a =2 (P), l($£; J=2 03Mi}| 2 , and we found that the 
maximum value of the squared overlap is approximately 
90 %. 

In the cases of the third and fourth states, ipm = 
*A=K^n) with tf ln =0.00165, and Viv = *£ 2 (^=0), as 
we mentioned already, there exist mixings of continuum 
states into these states. FIG. 7(b) tells us that the con- 
tinuum states which mix into ipiy seem to be dominantly 
the fourth lowest continuum state which we denote as 
ifi(6=0). But there also may be small mixing of the 
third and fifth continuum states, <f 3(6=0) and ip$ (6=0), 
into ipiv- In the case of the third state, ipm, there may be 
small mixing of the second and third continuum states, 
^2 (^111) and ^3(6111). Thus ipm and ipxv can be expressed 
as 
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FIG. 10: (a): Contour map of the squared over- 
lap, ($!^ J=2 (/3)|#^|(<5=0.00365)}| 2 , in the two-parameter 
space, /3x(= fly) and fl z . (b): Contour map of the squared 
overlap, |($3 a J=2 (/3)|*^z|(5=0)}| 2 , in the two-parameter 
space, /3 X (= fly) and fl z . Volkov No. 2 force is adopted. 



Here, </>r(<5iii) and cp R (6=0) are just the resonance states 
which we want to extract from ipm and ipw, respectively. 

The extraction of these resonance states, 4>r(6ui) and 
<Pr(6=0), can be made in the following way. We explain 
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the procedure in the case of <^r(<5=0). First we note 
that the three eigen-energy states of GCM, ^^z|(5=0), 
tyJ= 2 (S=0), and ^fz 2 (S=0), which are the 5=0 states on 
the third, fourth, and fifth lowest continuum trajectories, 
respectively, and which are marked by crosses on the <5=0 
vertical dashed line in FIG. El can also be expressed by 
a linear combination of four states, 0r(5=O), <p4(8=0), 
ip 3 (S=Q), and (p 5 (6=0), just like ^=1^(5=0). This 
means that (/>r(<5=0) can be expressed by a linear combi- 
nation of four GCM eigen-states 



tion <^ a J - 2 (/3), namely |($^ J -^(/3)|0 R ( ( 5=O)| 2 . We see 
that the maximum value of the squared overlap is more 
than 88 %. We further calculated the squared overlap 
of <j>n(d=0) with PI =2 (/3q)^^ 2 which is the projection 
operator defined in Eq. 
more than 90 %. 



N,J=2, 



MS = o) 



7 



^D\^{=\8 

A=4 



0). 



(41) 



Since the rms radius of 0r(<5=O) should be much smaller 
than the rms radius of any of the three continuum states, 
f4(S=0), (fi3(5=0), and (p5(8=0), the coefficients {D™} 
of the above equation should satisfy that the expression 



A V *A =2 



(5 = 0) 
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A=4 



J2 D 

A=4 



A V *A =2 



(6 = 0) 



(42) 



takes the smallest value possible. This condition deter- 
mines the coefficients {D 1 ^} and hence </>r(<5=0). Sim- 
ilarly, we can extract the resonance state <^r(<5iii) by 
linearly combining the GCM eigen-states, ^^=3(^111), 
tkii=^i=l(8ui), and ^aSK^ii)- The calculated 
squared overlap, |(^r(<J=0)|^v>| 2 = |Ci V | 2 = |£>a=6I 2 and 

I <<^R (fill) I VlII> 1 2 = j Cl 11 1 2 = I ^A=4 1 2 is °- 811 and °- 967 ' re " 

spectively. This means that the degrees of mixing of con- 
tinuum states into iprv and ipm are around 19 % and 3 
%, respectively. We have checked by using this technique 
that the two states, ipn and ipi, are never influenced by 
their neighboring crossing points. 

In TABLE. II I II we show various properties of the 
four resonance wave functions along the 2jj~ trajectory, 
tpi, tpn, 0r(5ih), and 4>r(S=0). We have checked that 
the energies of these four resonance states (namely the 
expectation values of the Hamiltonian H'(6) by these 
four resonance states) measured from the 3a threshold, 
E(8)—Efa(6), lie almost completely on the curve of the 



ACCC real energy for the state given in FIG. 8(a) 
This result shows the high reliability of our extracted res- 
onance wave functions of </>r(5iii), and </>r(6=0). We see 
that the rms radius of the resonance state becomes larger 
as the excitation energy measured from the 3a threshold 
becomes higher from ipi to </>r(5=0). 

The resonance state 0r(5=O) is our approximate wave 
function of the 2\ state. Since the excitation energy of 
this state measured from the 3a threshold is quite a bit 
higher than that of the 0^" state, its rms radius is much 
larger than that of the Oj" state unlike for the case of 
tpl. However, this state preserves the same character as 
the ipi state that has a very similar structure to the 0^" 
state of Bose-condensate character. It is shown in FIG. 



10(b) where there is displayed the squared overlap be- 
tween 0r(5=O) and the single 2 + condensate wave func- 



(|37|l . The maximum value is 
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FIG. 11: The energy eigenvalues for J w = 2 + of the 
Hill- Wheeler equation, Eq. I|25|l . ploted as a function of 
the coupling constant, 8. *(zf(<5=0), *^=|(<5=0), and 
$^7(5=0) are marked by crosses in order from the bottom, 
and *x=l( 5 in) and #^=5 (5m) by open circles. Volkov No. 2 
force is adopted. 
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FIG. 12: Contour map of the squared overlap, 
I ($3^ J=4 (/3) 1*3^2 (5=0.015)) | 2 , in the two-parameter space, 
Px(= P v ) and /3 Z . Volkov No. 2 force is adopted. 



The fact that the r.m.s radius of the calculated 2^ 
state has the extraordinary large value of 5.40 fm deserves 
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TABLE III: Expectation values of kinetic energy, nucleon-nucleon interaction energy, and Coulomb interaction energy of the 
relative motions between a particles, denoted as T — Tq, Vn, and Vc, respectively. Internal contributions of the a particle are 
subtracted. Calculation is made for the state by using Volkov force No. 2 and various values of the coupling constant 5. 
When 5=0, 0.00165, 0.00265, 0.00542, the wave functions, (f> R (5 = 0), <fe.(<5m), *a=s(<5 = 0.00265), and ^{zl{5 = 0.00542) are 
adopted, respectively. Units are in MeV except for Ri. m . s whose unit is in mi. 
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special attention. It means that this state has about ten 
times the volume of the 12 C ground state. One may, 
therefore, consider it as an "alpha halo state". As we 
mentioned in subsection III Dl we will show in our next 
paper that our 2^ wave function is very similar to that 
of Ref. [j| by Kamimura. 

Finally we study the state obtained in the previous 
subsection. We estimate the 4^ wave function by using 
the energy eigen functions on the 4^~ trajectory shown in 
FIG. |7(c]) as was done for the 2^ state. FIG. H2l shows 
the squared overlap, | ($^ J=4 (/3)|^(z %(6 = 0.015)) | 2 . In 
spite of the fact that FIG. 7(c) shows that the energy 



(A) and (B), respectively. 



eigen function, ^x=i(^ = 0.015) is sufficiently far from 
the nearest crossing point with the continuum state, the 
calculated squared overlap never has any non-negligible 
amplitudes in the inner interaction region. The root 
mean square radius of this state is 19.9 fm, and this state 
can hardly be regarded as a physical resonance state. If 
we make the coupling constant S smaller than this value 
of 0.015, going towards S — 0, the wave function on this 
4 2 trajectory will have a larger radius than this unphys- 
ical value of 19.9 fm. That is why we hesitate to sup- 
pose that the obtained energy position and width of the 
4:2 state reasonably correspond to a physical situation, 
although the obtained energy position for the effective 
nucleon force, Volkov No. 1, is comparable with the one 
by Uegaki et al. 




0.02 



FIG. 13: Analytically continued complex energy functions for 
the ^" and 0^" states are drawn by superposing them to FIG. 
7(a)| (A) is the analytically continued complex energy func- 
state, and (B) for the 0^" state. Note that 



tion for the Ot 
k 2 

[6,6] 



and for^ yi 



for the O3" state and k 2 6 6 i 



and k 2 7 7 i for the 



0J state are converged, so that the two trajectories in each 
case of the 0^" and 0^" states are hardly distinguishable in this 
figure. Volkov No. 2 force is adopted. 



D. The Oj and 0+ states 



As we see in FIGs. 7(a)|7(c) shown in the previous sec- 



tion, IIII Bl it turns out that there exist several resonance 
states going under the threshold as S increases. For the 
+ and 2 + states we see three resonance states and for 
the 4+ state we see two, up to 6 = 0.02. For the + state, 
the 0^ and 0^" trajectories seem to cross mutually in the 
region with S smaller than 0.01. In FIG. El we show 
the complex energy functions analytically continued for 
both Ojj" and 0^ states in the case of Volkov No. 2 force. 
Throughout this subsection we will refer to analytically 
continued energy functions of the 0$ and 0^ states as 



We see that the trajectories of the real parts of (A) 
and (B) cross each other around 6 = 0.002 before reach- 
ing (5 = 0. And we also see that the trajectory of the 
real part of (A) is not straight but quite distorted in 
the region, 0.003 < <5 < 0.006. This unexpected situa- 
tion of the crossing of two lines might cause the distor- 
tion for the trajectory of the real part of (A). For (A), 
Re{kf 7J] (S = 0)) = 4.8 MeV and T = -2Im(k? 7J] (S = 
0)) ='9.2 MeV, and for (B), Re(fc| 8] (5 = 0))'= 1.9 
MeV and V = -2Im(fc| g] (5 = 0)) = 0.36 MeV. This 
kind of situation including the non-straight shape of the 
real part of (A) and the crossing of the real parts of (A) 
and (B) similarly occurs in the case of the Volkov force 
No. 1. In this case, the complex energy function corre- 
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sponding to (A) gives Re(fc| g j(<5 = 0)) = 6.9 MeV and 
T = -2Im(fc| 9] (<5 = 0)) = 6.2 MeV. Both kf 7a] (S) and 
k^ s 8 i (5) give almost the same trajectory as fcp g g j (S) in a 
sufficiently wide energy region. And the complex energy 
function corresponding to (B) gives Re(£^ 8 g j(<5 = 0)) = 
1.8 MeV and T = -2Im(fc| g] (S = 0)) = 0.4 MeV. Again 
both k'^ 6 6 j (8) and k? 7 7 i (<5) give almost the same trajec- 
tory as fcp g g j (<5) in a sufficiently wide energy region. In 
order to determine the coefficients of all these Pade ap- 
proximants, we adopted the energy region from several 
tens of keV to around 3 MeV below the 3a threshold. 
Since the observed binding energy measured from the 3a 
threshold and total decay width of the Of state are 3.0 
MeV and 3.0±0.7 MeV, respectively, neither (A) nor (B) 
reasonably explains the experimental data. 

We wonder if the ACCC method where the analytic 
continuation is made using the Pade approximant works 
well in the case where this kind of crossing occurs. Also 
the distortion of the trajectory of the real part of (A) 
makes us feel doubtful whether (A) describes the reso- 
nance reasonably well. 

In section^ we referred to Refs. ^lj an d El m which 
the density distribution of the intrinsic state of the Of 
state is reported to have linear-chain-like configuration 
of 3a's. If this result is reasonable, our present approach 
in which strongly prolate deformations with (3 x =(3 y <0.b 
fm and (3 Z >3.0 fm are not included in the GCM calcu- 
lation may not be suited for the description of the Of 
state. However it is not clear at all whether this possi- 
ble inappropriateness of our approach is related or not to 
the crossing problem of the ACCC trajectories or the dis- 
torted shape of the ACCC trajectory of energy functions. 
Further studies are necessary to clear up the situation. 



IV. DISCUSSION 



Angular momentum component in the intrinsic 
wave function 



Almost all of the studies of this paper are based on 
the wave function with good angular momentum quan- 
tum number, Eq. (|13f) . which is projected out of the 
deformed intrinsic wave function, Eq. JT2). The de- 
formed wave function, $3 Q (/3), is characterized as form- 
ing a gas-like 3a-cluster structure. However, in order 
to be able to impose the same character as the intrinsic 
wave function on the angular-momentum-projected wave 
function, it is necessary that the projected wave function 
is contained in the intrinsic wave function with a non- 
negligible amplitude. We define ^ a (f3) as a normalized 



wave function of Eq. (|12fl . For J v = + and 2 + , we cal- 
culated the magnitude of the squared overlap amplitude, 
|($^ J (/3)|$^ (/3))| 2 , where (0) is given in Eq. 
The c ontour maps are shown in the FIGs. 14(a)| and 
14(b) for J 71 " = + and 2 + , respectively. The harmonic 



oscillator size parameter, b = 1.35 fm is adopted. Even 
for the 2 + state, a non-negligible amplitude can be seen 
in both prolate and oblate deformed regions. 



B. Correlation between the Of and 0^ states 

As already mentioned, the fact that the Of state has 
a dramatically different structure from the shell-model- 
like structure of the Of state was revealed by the works 
made about a quarter century ago [J 0- Q There, the 
cluster structure of the Of state was characterized by the 
feature that the wave function has a large amplitude in 
the outer region, while the amplitude in the inner region 
is suppressed. The formation of the clustering structure 
of the f state was attributed to the orthogonality to the 
Of wave function that has a compact structure, therefore, 
a large amplitude in the inner region. If this interpreta- 
tion of the appearance of the cluster structure of the Of 
state is correct, the structure of the Of state would not 
be affected by the detailed structure of the Of state if 
only the 0^~ state has a compact structure. 

In this subsection, we show quantitatively that the 
detailed structure of the Of state does not affect the 
structure of the Of state that is orthogonal to the Of 
state. In FIG. 1 of Ref. Q, we have shown the contour 
map of the energy surface of the + state in the two- 
parameter space, P x (= (3 y ) and (3 Z , obtained by adopt- 
ing the Volkov No. 2 force, which showed the minimum 
energy state at (3 X (— [3 y ) = 1.5 fm, fi z = 1.5 fm, i.e. 
/3o = (1-5 fm, 1.5 fm). The obtained minimum energy 
was —87.68 MeV, which was considered to correspond 
to the approximate binding energy of the Of state. In 
the reference paper, we further introduced the projec- 
tion operator, Pj 7=0 (/3o) which is shown in Eq. I|35() 
of this paper, and then looked for the minimum en- 
ergy state for the + state in the orthogonalized space, 
Pj 7=0 (/3 ) < i ) 3 < ^ (/3). The minimum energy state was ob- 
tained at Px(= Py) = 5.7 fm and (3 Z = 1.3 fm, and the 
minimum energy was —81.55 MeV, which was considered 
to correspond to the approximate binding energy of the 
Of state. 

In order to study the sensitivity of the Of state on 
the 0^~ structure quantitatively, we introduce /3o' = 
(1.0 fm, 2.0 fm) and calculate the energy surface us- 
ing the orthogonalized wave function Pf =0 (flo')®^ (fl) ■ 
The binding energy of the wave function, $3= (/3 ') is 
calculated as —85.96 MeV, using the Volkov force, No. 
2, which is higher by 1.72 MeV than the value corre- 
sponding to $3=°(/3 )- 

FIG. El is the contour map of the energy surface 
corresponding to the + state, in the orthogonalized 

space, Pj 7=0 (/3o')*3r (/ 3 )- We can see tliat tliere 
appears the minimum energy point at [3 x —[3 y —d>.2 fm, 
f3 z = \A fm with the minimum energy —81.41 MeV and 
the qualitative feature is quite similar to FIG. 2 of Ref. 
• In spite of the fact that the binding energy given by 
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(a) 



J K =2 + components 




(b) 



FIG. 14: (a): Contour map of the squared amplitude, 

K$3a J= V)l$3a(/3))| a in the two-parameter space, X (= f3 y ) 
and j3 z . (b): Contour map of the squared amplitude, 

|($3a J=2 (/3)|*8a(/3))| 2 in the two-parameter space, /3 X {= /3 y ) 
and j3 z . Harmonic oscillator size parameter 6=1.35 fm is used 
in both figures. 




FIG. 15: Contour map of the energy surface corresponding to 
the + state orthogonalized to the state at f3 o '=((3 x =/3 y =1.0 
fm, /3 Z =2.0 fm). The adopted effective nucleon force is Volkov 
No. 2. Numbers attached to the contour lines are binding 
energy values given in unit of MeV. 

<&3~°(/3o') is higher than the value given by c E > 3 t T°(/3o) 
by 1.72 MeV, the minimum energy and the minimum 
position in FIG. 1151 are almost unchanged compared 
with those calculated using Pf =0 (/3 )$3= °(/3). This 
result is consistent with the argument qualitatively made 
by the previous works referred to above, and further 
assures that the Oj" state could be separated from the 
§i state and treated independently under the condition 
that the state cannot become compact like the state. 



C. Dominance of spherical component in the 
condensate wave function 

If our GCM wave function of the 0^" state, ^x=2 ^ s 
extended to a deformed type of the a condensate, one 
may ask the question how important this effect is with 
respect to the spherical case ? 

The discussion of the III CI subsection teaches us that 
the deformed condensate wave function, after angu- 
lar momentum projection, necessarily contains a large 
amount of the spherical condensate component. Actu- 
ally we showed that after the projection onto the J T = + 
state, both prolate and oblate condensate wave functions 
take the form of Eq. (|19fl for the leading term which is 
nothing but the spherical condensate wave function. 

In order to check the above argument quantitatively, 
we investigated the magnitude of the spherical conden- 
sate component which incorporates ^fz®- Let ^sph an d 
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Wdef be the functional spaces spanned by spherical and 
by deformed condensate wave functions, respectively; 
V sp h Q Wdef, an d "'"Vsph be the orthogonal functional 
space to V sp h- The projection operator P sp h onto V sp h is 
written as P sph = £„ l*? P h)(*" P hl» wh ere *^ ph are the 
orthonormal basis functions of the space V sp h- The ^s ph 
are constructed as follows, 

£fe°(/% = $y = =Py = 

0a, 



for the + state, and 



sph 



^X> W (&)*£ (&=A, = &)• 



^a=2 can be expanded by a linear combination of 
*" p f, *" p l 2 ■ • • € V sph and a vector x € ^Vsph H W dc f as 



■ sph 



'sph "1 r X, 

where the relation, P sp hX — is satisfied. The 

spherical condensate component can be calculated as 
(*^j=0|Psph|*A=2 = °)- When thc number of the adopted 

components *™ ph in P sp h is 16, where (3 X = (3 y = (3 Z 
varies from 0.5 fm to 15.5 fm by 1 fm steps, the ob- 
tained values of the spherical condensate component 
are 92.1 % and 90.6 % for the Volkov force, No. 1 
and No. 2, respectively. Even when we decreased the 
number of the adopted components to 13 from 16, i.e. 
Ac = fly = f3 z = 0.5 fm to 12.5 fm by 1 fm steps, there 
appeared no difference of the obtained values at least to 
six significant figures for both forces. From this result it 
follows that our obtained values have already converged. 
The large magnitudes of these values imply that 
is composed of the spherical condensate component by 
more than 90 %. At the same time we have to note that 
some amount (less than 10 %) of the deformed compo- 
nent orthogonal to the spherical one is necessary in order 
to have quantitatively good reproduction of the observed 
results. 



D. Vibrational structure of the angular momentum 
projected wave function of the deformed condensate 

The expressions of the leading terms of the angular mo- 
mentum projected wave functions of the deformed con- 
densate, Eq. ljl§l) for the 0+ state and Eq. lj2*U|) for the 
2 + state, can be rewritten as follows, 



n-l 



A 



exp {- } </>"(") 



= A 



9 2 



(43) 
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n 

x J2( X i - X G ) 2 Y 2Q (Xi^X Gl W{ a ) 



i = l 



n 

= Y t A[{X j - X G ) 2 Y 20 (Xi^X Gi ) 

3=1 

n 2 
xI1 cx p{- z*(Xi-XG) 2 }<t>(<*i)], (44) 
i=i 

for the 2 + state. The functional form of Eq. (|43|l means, 
as was mentioned several times, that all n a-particlcs 
occupy an identical S'-orbit, exp{ — (2/B 2 )(X — Xq) 2 }. 
On the other hand, the functional form of Eq. I|44l) 
means that one a-particle occupies a -D-orbit, (X — 
X G ) 2 Y 20 (X - X G ) exp{-(2/P 2 )(X - X G ) 2 } while the 
other (n — 1) a-particles occupy an identical S-orbit, 
exp{-(2/P 2 )(X-X G ) 2 }. Therefore, our 2+ wave func- 
tion constructed by the angular momentum projection of 
the deformed condensate wave function has the structure 
that one a-particle is excited from the S'-orbit into a D- 
orbit while the other (n—1) a-particles remain in the 
identical S'-orbit which is occupied by all n a-particles in 
the case of the + state. 

We can yet give another explanation for the 2 + wave 
function and this explanation can be applied to any wave 
function with spin obtained by the angular momentum 
projection of the deformed condensate wave function. We 
rewrite Eq. (|18fl as follows, 



ll—x 

exp | - 2 m ( 



t 2 +£ 2 £ 2 



m 



B* 



)} 



i=l 



exp 
xjl + oD 



° 2 ^ , - 3 



-rD J + ^— D 3 
2! 3! 



}, (45) 



D = ^(X 4 - X G ) 2 Y 2Q (X t ^~X G ), 

i=i 

2 fldn / 1 1 \ 

We can say that the leading term of the 2 + state is cre- 
ated from the + state by one P/-phonon-like excitation 
expressed by the operator D. Similarly, the leading term 
of the 4 + state is created by two P-phonon-like excita- 
tion expressed by the operator (P/ 2 )j = 4, where (P 2 )j = 4 
stands for the operator obtained from D 2 by the angu- 
lar momentum projection onto J — A. We can easily 
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see that any spin J wave function obtained from the de- 
formed condensate wave function has the leading term 
which has the form of multiple Z?-phonon-like excitation 
from the + state expressed by the operator (D m )j. 

The above argument shows that all the angular mo- 
mentum projected wave functions of the deformed con- 
densate can be viewed to express the state of vibration- 
like excitation rather than the rotational excitation. 
However, the GCM calculation by superposing the spin- 
projected condensate wave function is not directly re- 
lated to the diagonalization of the Hamiltonian whithin 
the multi-phonon-like subspace with definite J, because 
in the GCM calculation we superpose phonon-like states 
with different oscillator parameters B. 



V. CONCLUSION 

We investigated the states with J 7T =Q + , 2 + and 4 + , of 
12 C with the excitation energy less than about 15 MeV 
by using the 3a condensate wave function with spatial 
deformation. In order to discuss the total alpha decay 
widths as well as the energy positions we applied the 
ACCC method to the wave function. 

We could successfully reproduce both the binding en- 
ergy and total alpha decay width of the recently observed 
2\ state. The observed values of the energy and width 
are 2.6 ± 0.3 MeV and 1.0 ± 0.3 MeV, respectively. The 
calculated binding energy and width are 2.1 MeV and 
0.64 MeV, respectively, for the use of Volkov No. 1 force. 
For the use of Volkov No. 2 force, when the binding en- 
ergy is adjusted to 2.0 MeV with the nuclear interaction 
being a little weakened, the width is calculated as 0.61 
MeV. We found that in the case where the binding en- 
ergy of the 2\ state is adjusted to 0.26 MeV which is the 
calculated binding energy of the 0^", the density distri- 
bution of the 2^~ state is almost the same as that of the 
0^~ state for the use of Volkov No. 2 force. And then, the 
kinetic energy, nuclear interaction energy, and Coulomb 
interaction energy of the calculated 2^ state are found to 
be very similar to those of the 0^" state. 

We devised a new technique to extract the proper reso- 
nance component from the wave function with some mix- 
ture of continuum states which are obtained within the 
bound state approximation. This technique allows us to 
obtain pure resonance wave functions very easily within 
the framework of the bound state approximation. De- 
tailed discussions about this new method will be given 
in our forthcoming paper. As a result of applying this 
method to our present problem, we found that the pure 
resonance 2 J wave function which gives almost the same 
resonance energy as the above mentioned ACCC result 
has a large overlap with the single condensate wave func- 
tion of 3a gas-like structure. The squared overlap value 
is more than 88 % when we employ the Volkov No. 2 
force. 

These facts imply that the 2 J state has a similar struc- 
ture to that of the 0^ state which has a gas-like struc- 



ture composed of three alpha clusters and is of Bose- 
condensate character. The 2^" state is obtained by pro- 
moting just one alpha cluster out of the condensate of the 
O2 state into a D-w&ve. The result of the present paper 
about the 2^~ state of 12 C suggests that for more heav- 
ier nuclei, for instance, 16 O, there would exist non-zero 
spin excited states which are members of an alpha cluster 
condensate family. The 0^" state is a state of very dilute 
density which is about 1/3 of the ground state density. 
The density of the 2% state is pointed out to be even more 
dilute than that of the 0^" state. The 2 J state has an ex- 
traordinarily large r.m.s radius such that the volume of 
the l\ state is a factor ten larger than the one of the 
ground state. We call such a wide a-configuration an a- 
halo state. These dilute densities or large radii of the 0^ 
and 2^~ states are probably the reason why no-core shell 
model calculations of 12 C completely fail to account for 
the O2 and 2;j~ states in spite of their good reproduction 
of the ground state band [24[ . All these above features 
hint to the possibility of a description of these states in 
terms of ideal bosons as this has in deed been considered 
in @. 

We failed to give a clear interpretation of the 0|" state. 
As mentioned in the text, two recent works [l^, E3 ob- 
tained a practically identical structure for the 0jj~ state 
consisting out of a 8 Be-like entity with a third a attached 
to it by about 30° off the symmetry axis of the first two 
a's. Though these calculations fail to reproduce the 0^" 
position by several MeV, such results are nevertheless of- 
ten indicative of the geometrical structure of the state 
as this is also the case for the 0^" state in 12 C obtained 
from an AMD aproach In these calculations, the 0^" 
state looks similar to a chain state but not quite. May 
be a very prolate deformed a-condensate wave function 
could elucidate the physical interpretation of the state. 
It is our next task to investigate the 0^ state by incorpo- 
rating strongly prolate condensed wave functions in the 
GCM calculation. 
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APPENDIX: MATRIX ELEMELNTS OF 
DENSITY OPERATOR 

We explain how to calculate the matrix elements of the 
density operator defined in Eq. Since the density 

operator maintains the rotational invariance, the expec- 
tation value by the normalized wave function with good 
angular momentum can be written as follows; 





p(a) 




(S 


^na 


3))( 
p(a) 




) 









(A.l) 

In order to evaluate the above expression, we make use 
of a similar technique to the one used in Eq. I|15[l. At 
the final step of Eq. I|15|) . we change the numerator as 
follows, keeping the denominator unchanged: 

d cos 6d J Q0 (9) J^^nc ((3) | p(a)R y (9) | <f„ Q ((3'] 



= d 



(*n«G9) 


S(X G )p(a)R y (9) 


$n«(/3')) 


(«f£ G (/3) 


S(X G )R y (9) 


<fr£ G (/3')) 



dcos0<tf o (0)(* na (/3) 6(X G )p(a)R y (9) $„ Q (/3') 



Here, we made use of the relation of Eq. Ullfl 

$„„(/3) = Q$* G (/3) 8„ a (/3), 
where Q is a constant, $^ G (/3) is defined as 

dr^(/3) = exp(- £ 

k—x,y,z 

and use is made of the trivial relation 



(/3') 



$^ G (/3) 6(X G )Ry(0) 
Next, by using the following relation, 

J d 2 aY / Y l * m (r t ^X G )Y lm (a) = 1, 



(A.2) 
(A.3) 

(A.4) 

(A.5) 
(A.6) 



I III 



the density operator is represented below as 

i A 

P(a) = - YS(\ ri -X G \-a) 



A 



2 A r 

^jd 2 a%-J G -a). (A.7) 



Here a and r, — X G are polar angles of a and — Jg, 
respectively. Substituting Eq. (|A. T|) into the final expres- 
sion of Eq. (|A.2|I , we only have to evaluate the following 
expression, while the integral over the variables, cos 9 in 
Eq. l|A~2|) and a in Eq. (|A.7|I . is numerically performed 
afterwards. 
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Considering the fact that the wave function, $ na (/3) is 
represented by using Brink's wave function, as 
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we should be able to evaluate the following quantity, 

A 

Y8(ri-a)\$ B (R' 1 ,---,K) 

i=l 

d 3 k 



(2tt) 3 



$ B (Hi, •••,#„) JJ exp(-fc-r J 



5>(n-a) $ b ( j r' 1 ,--., j r: i 

1=1 

(A.10) 

All of the following procedures are well known and can be 
straightforwardly performed, once a single particle wave 
function, ifi(r), and a modified single particle wave func- 
tion, ip'^r), are introduced as follows, 

Vi (r) = (^ 2 )- 3 / 4 exp(-^(r- J R,) 2 ), 

^(r) = (^ 2 )^/ 4 exp(-^(r- J R' t ) 2 + ^fe.r), 

(i = 1, ■■■,n). 

(All) 
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The overlap kernel between ipi(r) and f'Ar) is calculated 



as 



(A.14) 



Di 
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(A.12) 



where we use a shorthand notation, = (ipi\ipj). Z)y 
and are defined as n x n matrices, and the inverse 
matrix of Z?y is also calculated easily, 

= ex P - ^4 fe ' & + 

(A.13) 



Now we will see that using Eqs. (|A. 12|> . (|A.13|I . Eq. 
(|A.10(1 is represented as follows: 



f d 3 k I | JL » 

/ (^r^'-'-'^ln ex P (-fe-r,)x 

i=i 



(2tt) s 



4|D| 4 £(^|5(r--a)|^)(5- 1 ),- i . 



After substituting the explicit formulae of Dij and 
(D~ 1 )ji into the above expression and performing the 
Gauss integral over the variable, k, analytically, we fi- 
nally arrive at the following formula: 



exp 



Ri + b!a \ 2 



6 2 (A- 1) 



r v— -\ , / , R% + R'i ~i 2 

-{a + 2Y / (Rk + R' k ) 2~^} 



(A.15) 



Note that it is not complicated to perform the Gaussian 
integral over the variables, R\, ■ ■ ■ , R n . As mentioned 
above, the integral over the variables, a and cos#, is per- 
formed numerically. 
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